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ABSTRACT 

An unsolved problem in turbulent dynamo theory is the "back reaction" problem: to 
what degree does the mean magnetic field suppress the turbulent dynamo coefficients which 
are needed to drive its growth? The answer will ultimately derive from a combination of 
numerical and analytical studies. Here we show that analytic approaches to the dynamo and 
back reaction problems require one to separate turbulent quantities into two components: 
those influenced by the mean field (which are therefore anisotropic) and those independent 
of the mean field (and are therefore isotropic), no matter how weak the mean field is. Upon 
revising the standard formalism to meet this requirement, we find that: (1) The two types of 
components often appear in the same equation, so that standard treatments, which do not 
distinguish between them, are ambiguous. (2) The usual first-order smoothing approxima- 
tion that is necessary to make progress in the standard treatment is unnecessary when the 
distinction is made. (3) In contrast to previous suggestions, the correction to the dynamo 
a coefficient found by Pouquet et al (1976) is actually independent of the mean field, and 
therefore cannot be interpreted as a quenching. 

Subject Headings: magnetic fields; galaxies: magnetic fields; Sun: magnetic fields; 
stars: magnetic fields; turbulence; accretion discs. 



- 2 - 



1. Introduction 

Mean-field turbulent dynamo theory is a leading candidate to explain the origin of 
large-scale magnetic fields in stars and galaxies (e.g. Moffatt 1978; Parker 1979; Krause 
& Radler 1980; Zeldovich et al. 1983). This theory appeals to a combination of helical 
turbulence, differential rotation, and turbulent diffusion to exponentiate an initial seed large- 
scale magnetic field. Typically, the total magnetic field is broken into a large-scale or mean 
component and a small-scale or fluctuating component, and the rate of growth of the mean 
field is sought. The mean field grows on a length scale much larger than the outer scale of 
the turbulent velocity, with a growth time much larger than the eddy turnover time at the 
outer scale. Fluid helicity provides a statistical correlation of small scale loops favorable to 
exponential field growth. Turbulent diffusion is needed to redistribute the amplified mean 
field since Ohmic diffusion is negligible. By definition, the mean-field dynamo then belongs 
to the "slow" dynamo class (Zeldovich et al 1983). On the other hand, the small-scale 
field amplifies quickly without the aid of diffusion, and is thus constitutes a "fast" dynamo 
(Zeldovich et al 1983). Rapid growth of the small-scale field necessarily accompanies the 
mean-field dynamo. Its impact upon the growth of the mean field, and the impact of the 
mean field itself on its growth are controversial. 

Dynamo theory has usually been studied in the kinematic regime, wherein the turbulent 
and ordered velocity fields are assumed to be given, and the seed field is assumed to grow 
and amplify without back reactions that affect the turbulent motions which drive its growth. 
Kinematic dynamo theory is incomplete, as a real dynamo should represent a solution to the 
nonlinear magneto-hydrodynamic (MHD) equations, which certainly embody back reaction. 
Although numerical simulations are necessary for our understanding of the nonlinear theory, 
the simplicity of the basic mean-field dynamo formalism warrants analytic investigation of 
the nonlinear MHD equations and the resulting back reaction to understand why or why 
not something like the kinematic treatment captures the evolution of the field even in the 
nonlinear regime. 

Some aspects of how the growing magnetic field affects the turbulence have been ad- 
dressed. Kulsrud & Anderson (1992) confirmed previous arguments by, e.g., Cowling (1957) 
& Piddington (1981) that the small-scale field builds up to approximate equipartition suc- 
cessively from the smallest to the largest scales in the turbulence. However, since most of 
the energy in a turbulent spectrum is concentrated at the largest - or outer - scale, the 
question of whether the outer scale reaches true equipartition is most important, since at 
equipartition, MHD turbulence consists largely of Alfven waves which do not contribute to 
dynamo action (Moffatt 1978); this would strongly suppress dynamo growth. An equally 
important, and related question is: for what value of the mean field does suppression of dy- 
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namo action become important? The simulations of e.g. Cattaneo & Hughes (1996) suggest 
that suppression occurs when the energy density in the mean magnetic field is smaller by a 
factor of the magnetic Reynolds number than that implied by e.g., Pouquet et al. (1976). 

Not only do simulations disagree on the role of the magnetic Reynolds number in the 
suppression formulas, but so do interpretations from well-motivated analytical approaches 
such as those of Zeldovich et al. (1983); Montgomery & Chen (1984); Gruzinov & Diamond 
(1994); Bhattacharjee & Yuan (1995); Kleeorin et al. (1995) and Field et al. (1998). The 
formalism of Bhattacharjee & Yuan (1995), Gruzinov & Diamond (1994) and Kleeorin et 
al (1995) leads to extreme suppression even in the weak mean-field limit. A second class 
of formalism is exemplified by Vainshtein & Kitchatinov (1983) and Montgomery & Chen 
(1994) who considered an expansion in the mean magnetic field, and Blackman & Chou 
(1997) who considered a similar perturbation theory in both the mean magnetic field and 
the mean velocity field. Field et al. (1998) developed a fully nonlinear approach for assessing 
the back reaction in the case that the gradient in the mean field is negligible. In contrast to 
Bhattacharjee & Yuan (1995), Gruzinov & Diamond (1994) and Kleeorin et al. (1995), but in 
agreement with Kraichnan (1979), Field et al.(1998) find that the magnetic Reynolds number 
does not strongly enter the suppression of the alpha dynamo coefficient. Why do different 
analytic approaches in such papers concerned with back reaction produce such different 
results? Here we show that the answer to this question depends on how one decomposes 
small-scale turbulent quantities into components which are dependent, and independent 
of the mean fields, respectively. We show that these components often appear together 
in the same equations and not distinguishing them leads to ambiguities, while properly 
distinguishing between them allows progress even if the first-order smoothing approximation 
is not satisfied. 

In section 2 we review the standard derivation of the dynamo coefficients following Mof- 
fatt (1978). In Section 3 we dissect the standard treatment, point out the ambiguities, and 
rederive the a coefficient with a formalism that makes the distinction between components 
discussed above. In section 4 we show the effect of the revised approach on analytic back 
reaction studies. We conclude in section 5. 



where A is a constant magnetic viscosity. Here B is in velocity units, obtained by dividing 
by \/4:irp; we assume incompressibility. 



2. The standard derivation 



The induction equation describing the magnetic field evolution is 




(1) 
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The equation for the mean field derived by averaging ([!]) is 

d t B = V x (v x b) - V • VB + B • VV + AV 2 B, (2) 

where an overbar indicates the mean field and the small-scale quantities are in lower case. 
The term B • VV describes the effects of differential rotation, and will not be discussed 
further here, while the term V • VB can be eliminated by changing the frame of reference to 
one moving with V; both terms will be ignored in what follows. The task for the dynamo 
theorist is to find the dependence of the turbulent EMF (v x b) on B so that (§) can 
be solved. The kinematic theory assumes that the statistical properties of the small-scale 
velocity field are prescribed independently of the mean magnetic field; in that case it remains 
only to obtain the statistical behavior of the fluctuating magnetic field. Subtracting (0) from 
(0) then yields 

d t h = B ■ Vv - v • VB + b ■ Vv - v • Vb - V x (v x b) + AV 2 b. (3) 

At this stage, it is customary to make the first order smoothing approximation (FOSA), 
namely, that the terms second order in the fluctuating quantities can be neglected. The 
result is then 

d t b = B • Vv - v • VB + AV 2 b, (4) 

and for small A the last term can be ignored in the kinematic theory. We use this equation 
to describe the statistics of b in terms of those of v. In doing so we employ the Reynolds 
relations ( |Radler (1980)|) , i.e. that derivatives with respect to x or t obey d t ^(XiXj) = 
(df^XiXj)) and (XiXj) = where X{ = Xj + X, are components of vector functions of x 
and t. For statistical ensemble means, these hold when correlation times are small compared 
to the times over which mean quantities vary. For spatial means, defined by (Xj(x, i)) = 
V^ 1 J Xj(x + s, t)ds, the relations hold when the average is over a large enough V that 
L < 1/5 <C D, where D is the size of the system and L is the outer scale of the turbulence. 

Using the Reynolds rules, isotropy of the turbulent velocities, and the assumption that 
(v(t) x b(0)) = 0, we obtain the classical result (Moffat 1978): 

(v x b) = (v(t) x C d tl bdt') = aB — /3V x B, (5) 
Jo 

where a = — |t c (v- Vx v) and (3 = |t c (v-v) are the helicity and diffusion dynamo coefficients 
respectively, and t c is the correlation time of the turbulent quantities. Eq. @ is put into 
(0) to solve for the evolution of the mean field, often assuming Va = V/3 = 0. 

3. Including the back reaction and improving the standard formalism 
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There are three fundamental limitations with the derivation of section 2: (a) The kine- 
matic approximation, according to which the velocity field is prescribed, ignoring the back 
reaction of both the small- and large-scale magnetic field, (b) The FOSA, according to which 
1st and 2nd terms in (3) dominate the 3rd and 4th terms, even though it is known (Kulsrud 
and Anderson 1992) that the small-scale field b builds up much faster than the large-scale 
field B. (c) Finally, even though the same symbols v and b are employed throughout the 
derivation, they mean different things in different places. To see why, note that the scalar 
forms of a and /3 follow from isotropy of the fluctuating components b and v. However, if 
both are strictly isotropic, (v x b) would vanish because it is the average of an isotropic 
vector. Moreover, growth of the mean field in ([|) requires some violation of homogeneity, 
because growth requires a non-zero V x (v x b). Thus, some way of distinguishing the small- 
scale quantities in the turbulent EMF from those in terms of which it is expanded using (TJy 1 
is necessary. 

We will first formulate the required modifications to the standard theory. Then we 
study their effects in attempting to derive the Reynolds number dependence of the back 
reaction on a in the limit of weak mean field. To this end, we split up the equation for 
the small-scale magnetic field into an equation for that component which is independent of 
the mean magnetic field and that component which depends on the mean magnetic field. 
Blackman Sz Chou (1997) applied this approach to first order in both the mean field and the 
velocity field to find a coupled vorticity-magnetic dynamo. With a similar approach, Field 
et al. (1998) compute rigorously the effect of the back reaction of B on the dynamo a to all 
orders in B; here we consider only the first-order terms in their development. 

We assume that the turbulence is weakly anisotropic and inhomogeneous. Terms lin- 
ear in the slowly time-varying mean quantities contribute, but their averaged zeroth-order 
coefficients are taken to be isotropic and homogeneous (still allowing for reflection asymme- 
try). Iterating the equations using the formal solutions for the turbulent fields b(£) = h(t = 
0) + / dfhdt' and v(t) = v(i = 0) + / dfvdt', and using times appropriately chosen such that 
the correlation (v(£) x b(0)) ~ 0, we obtain to first order in mean quantities the turbulent 
EMF t f 

(v x b) (1) = (v (0) (t) x / <9 t ,b (1 W)-(b (0) (t) x / dt^dt'), (6) 

Jo ' Jo 

where the time derivatives are given by (3) and by the equation of motion, and thus to first 

order, (Q) becomes 

9(B = Vx(vx b) (1) + AV 2 B. (7) 

Here another ambiguity arises if one does not properly separate the zeroth-order quantities 
from the higher order quantities: Note that the quantity (v x b) can be written as (v x b) = 
(v(t) x / d t 'h(t')dt') = — (b(t) x J dt>v(t')dt'}. But not distinguishing between v, b and 
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v(°),b(°) in (|) leads to (v x b) = (v(t) x J d t b(t')dt') = ~{b{t) x / d t v{t')dt') = (v(t) x 
/ d t b(t')dt') — (b(t) x / dtv(f)) = 0, where the second last equality follows from purposely 
ignoring the superscript ^ in (6), and the last equality is therefore unavoidable. This shows 
that separating the zeroth order fluctuating quantities from the total small-scale quantities 
is essential. 

According to (6) the calculation of the EMF requires expressions for b^ and . The 
equation that determines is the momentum equation, given in general by 

d t V = -V • VV - Vp e ff + zA7 2 V + B VB + G + F(x, t), (8) 

Here F is the force that drives the turbulence, p e g = p + \B 2 , with p = P/p, v is a constant 
viscosity, and G is gravity. To put this equation in the form that we need, we first define 
the small-scale velocity v to be V — V and similarly for the other variables. Then averaging 
the resulting equation for V + v and subtracting the result from we obtain 

d t v = -V ■ Vv - v ■ VV - v ■ Vv + (v ■ Vv) 

-Vp- V|6 2 + V(\b 2 ) - V(b-B) (9) 
+b • Vb - (b ■ Vb) + b • VB + B • Vb + f (x, t) + iA7 2 v, 

where p is now the small-scale pressure and f is the small-scale (and only) part of F; consistent 
with the assumption of incompressibility, we take the small-scale gravity to be zero. We 
eliminate the term V • Vv by a change of frame. We shall also ignore the term v • VV. In 
our Galaxy for example, this term is smaller than |v- Vv| (Field et al 1998). In what follows, 
therefore, V will not appear. Now, writing the equations (|3]) and @ to zeroth and linear 
order in B, we have 

dth (o) = b (o) . Vy (o) _ y (o) . vb (o) _ v x ( v (o) x b (o)j + AV 2 b(°) (10) 

and 

d t bW = B • Vv(°) - v(°) • VB + bW ■ Vv(°) - v^) • VbW 
+b (°) . VvW - v(°) • Vb« - V x (v x b)^ + AV 2 b^). ( } 

for the small-scale magnetic field, and 

d t v^ = _ v (o) . v v (°) + (v(°) • Vv(°') - Vp (0) - V|6 (0)2 + b(°) • VbW - (b< > • Vb^) 

+ V(|6(°) 2 ) + f(x,t)(°) + z/V 2 v( ), 

(12) 

and 

^v^ 1 ) = -v(°) ■ VvW - v^) ■ Vv(°) - V^ 1 ) - V(b(°) ■ bW) + V(b(°) ■ b«) + (v • Vv)« - V(b(°) ■ B) 

+b (i) . V b(°) + b(°) ■ Vb« - (b • Vb)W + b(°) VB + B VbW + z/V 2 v( 1 ), 

(13) 
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for v 1 - 1 ) , where we have assumed that the forcing function is the same for both the state with 
B = and that with B ^ 0, so f (1) = 0. 

To further simplify we need an approximation related to, but distinct from the standard 
FOSA. Note that the standard FOSA, employed in section 1, would imply that |b| = |b(°) + 
b^ 1 ) + ...| < |B|. Although this assumption is different from the kinematic approximation in 
that it says nothing about the back reaction of the mean field on the velocity flows, it does 
imply that the small scale field is much weaker than the mean field. In reality, this is not 
the case: both simulations and analytical arguments (e.g. Parker 1979) show that the small- 
scale field energy saturates at a value of order the small-scale kinetic energy . Fortunately, 
we do not have to make the FOSA here. This is because we take the zeroth-order turbulent 
quantities to be solutions of (10) and (12), and work with the effect of the mean field on the 
first order quantities. We require only that |b^|, (v^l < B, a much weaker condition. This 
reduced smoothing approximation (RSA), unlike the FOSA, allows for » B. 

This RSA allows us to eliminate terms second order in the fluctuating quantities in 



(|T3|) and ([11]). Note, however, that in the first ~ t c of growth, the first term in ( pT3|) could 
dominate the terms with the mean field. However, we are only interested in times long 
compared to this time, after which a steady state ensues. When we put the equation for 
into the average (|6]) none of the bracketed terms on the right of (11) or (|T3|) contribute. The 
relevant equations for the small-scale fields then become 

= -Vp« - V(b<°> ■ B) + b(°) ■ VB + B VbW + zA7 2 v«, (14) 

and 

dtbW = B • VvW - v(°) ■ VB + AV 2 b«. (15) 

A remaining problem is to deal with the pressure in ([14]). For present purposes, we focus 
on the correction to the a dynamo coefficient, and so we consider the special case where VB 
is negligible (Field et al. 1998). We can then take the divergence of (14). The resulting 
equation is V 2 (p^ + b( ) • B) = 0, whose solution is (p^ 1 ' + b^ ^ ■ B) = constant. Thus the 
first two terms drop out of (14). 

Approximating time integrals in (v x b)^ by factors of the correlation time t c and freely 
employing Reynolds rules and incompressibility, we then obtain 

(v x b)W = a(°)B (16) 

where 



a 



(o) = -|t c [(v(°) ■ V x vW) - (b(°) ■ V x b(°))]. (17) 
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Upon substituting these into (fj), the curls can be pulled onto the B from homogeneity 
of the zeroth-order averages. Note that we have ignored the sub-dominant v and A terms. 
This will not pose any problem for our discussion in the next section: the entrance of the 
Reynolds number into the theory that we will address is that purported to enter through 
the "magnetic correction" (b^ ) ■ V x b^ )) in fll7|). Similar magnetic corrections to that in 
(|17|) have been derived before (e.g. Pouquet et al., 1976; Gruzinov & Diamond 1994) but 
none of the previous papers distinguish between a and aS Q \ and we will see below that the 
difference is fundamental. 

4. Contrast to previous work 

Building on the ideas of e.g. Piddington (1981), most workers agree that the Lorentz 
forces from the growing magnetic field react back on the turbulent motions driving the 
field growth and complicating the turbulent motions. However, analytic studies and sim- 
ulations disagree as to the extent to which the dynamo coefficients are suppressed by this 
back reaction. Some (e.g. |Cattaneo fc Vainshtein ( 1991 )| , [Vainshtein fc Cattaneo (1992 



Cattaneo ( 1994 )| , |Cattaneo fc Hughes (1996)] , Gruzinov fc Diamond (1994)| , Bhattacharjee 



fe Yuan (1995)| , [Kleeorin et al., (1995)|) argue that the suppression of e.g. a takes the 
form a ~ a^/(l + RmB 2 /v 2 ) where Rm is the magnetic Reynolds number, while others 
(e.g. Kraichnan 1979, Field et al., 1998) suggest a ~ a^°'/(l + B 2 /v 2 ) in the fully dynamic 
regime. If the former formula represented the actual level of suppression, the large values 
of Rm in nature would prevent astrophysical dynamos from working. Something else would 
be needed to generate large-scale fields in stars and galaxies. More analytic and numerical 
studies will ultimately be required to answer this question. At present, we suggest that some 
existing analytic arguments for extreme suppression can be challenged. 

To do so, we follow the path of pruzinov fc Diamond (1994)| and Bhattacharjee fc Yuan 



(1995)| but employ the formalism of section 3, rather than that of section 2, and then find 



that we do not arrive at the same conclusions. Consider the small-scale electric field: 

e = c 4 [(vxb)-vxb-vxB]+i)j, (18) 

where r] = AttX/c 2 is the resistivity. Dotting this with b and averaging gives 

( b . e ) = - C - 1 (b-(vxB))+r ? (b-j), (19) 

where we note that the first two terms on the right of flTS|) do not contribute to flT9|). Using 
the triple product on the second term in (^) and rearranging gives: 

(b-j) = r / - 1 [-B-(vxb/c) + (b-e)]. (20) 

Now the "magnetic correction" to a in the linear regime as computed in (17) is proportional 
to c(b- V x b)(°)/47r = (b-j)(°), not to (b-j). It is very important to distinguish between 
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those quantities with and without the superscript (0), as we now show. Taking the zeroth 
order contribution to fl2"0|) gives 



( b .j)(°)=r ? - 1 (b.e)(°), (21) 
but using fll8|) for e^ ^ in fl2~l|) shows that (^) is a trivial identity. Thus it is not possible 



to get information from (2~0p. If we do not distinguish the zeroth-order quantities from first- 
order quantities we might conclude that the left hand side of fl2~T|) is equal to the correction 
to a derived in (|TTD, but that is not correct. 



An almost identical discussion applies to Kleeorin et al (1995). In their Appendix A 
they provide a somewhat similar derivation of the back reaction to a, with the resulting 
equation 

d t a m = d(B ■ V x B - aB 2 /(3) - C 2 a m , (22) 

where C\ and C 2 are defined constants, and they write a = a +a m , with their a oc (v-Vx v) 
and a m oc (b ■ V x b). The subtlety arises from their statement below their un- numbered 
equation between their (A4) and (A5) that "a and t]t(= 0) are scalars when the mean 
field is small." Scalars result from the assumption of isotropy, but isotropy applies only 
to strictly zeroth-order quantities as discussed in section 3. Thus, the ambiguity in the 
notation discussed above is present. We suggest that the a and rj T {= (3) to which Kleeorin 



a 



et al (1995) refer should actually be and rfj> (= (3^)- Thus, the magnetic correction to 
in the weak mean field limit which they define as a m , should really be a$ oc (b^ • V x b^ )), 
not oc (b • V x b). This latter quantity is not simply related to a$ . The distinction is 
essential. Rather than (|22|), the equation for the correction would instead be 



ftog) = -Crf, (23) 

where no mean field appears. As a result, ( p2]) implicitly contains less information than it 
appears to, by analogy to (|20| - ^T1) . 



Yet another subtlety arises when one notices that the right side of ( pi]) is proportional to 
the sum of a time derivative and spatial divergence of averages of zeroth order correlations. 
Assuming exact stationarity and isotropy of zeroth order averages then implies (b • e)(°) = 0. 
The left side of (|2l|) would then also be zero, and thus the correction in ( plj ) would vanish 
for any finite r). The current helicity would then not be a correction to a but would be 
related directly through (pup. However, the fact that the diffusivity enters ( |2"T| ) means that 
the right side is the ratio of two small quantities and requires an assessment of how accurate 
the assumptions of stationarity and isotropy really are. 

The vanishing of (b • e)^ is reminiscent of Seehafer (1994) and Keinigs (1983), but is 
different because these authors do not distinguish zeroth order quantities and Keinigs (1983) 
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employs FOSA. Since averages of non-zeroth order fluctuating quantities are not in general 
isotropic, homogeneous or stationary, the full (b • e) ^ 0. When the full dynamo equation 
is considered, mean quantities that are functions of non-zeroth order fluctuating quantities 
can vary on the large scale of the system. The second term on the right of (E^) can exceed 
the left hand side, providing the dominant balance to the first term on the right. This is not 
apparent from Keinigs (1983) or Seehafer (1994) because they do not distinguish between 
zeroth-order and higher order small-scale quantities. 

5. Discussion and Conclusions 

We have pointed out a fundamental ambiguity in magnetic dynamo theory, resolved it, 
and have shown its importance in explaining why some analytic approaches to the problem 
of the back reaction of the mean field in dynamo theories disagree. In particular, we showed 
that only part of the small-scale quantities can be thought of as isotropic, no matter how 
small the mean field is, and why the formalism of standard treatments is ambiguous on 
this point. The "magnetic correction" term shown in (\Tf\j is constructed from only that 
component of the small-scale field which is independent of the mean field, not from the 
total small-scale field. Previous studies have not separated the small-scale quantities into 
their zeroth and higher order components, and as a a result find that the correction term 
in ( |TTD appears to relate to the magnetic Reynolds number through fl22|) or fl20|). When the 
separation is made, the correction term cannot be interpreted as a strong suppression to the 
a effect because then the only place for the magnetic Reynolds number to enter is through 
the last term in (|l"5|). This does not produce a strong suppression of a (Field et al 1998). 
Note that the first-order smoothing approximation is not necessary in our approach. 

The fact that a true dynamo cannot be kinematic need not imply that dynamo action is 
suppressed. An often quoted ratio, that of the total magnetic energy to that contained in the 
mean field (e.g. |Zeldovich et al. (1983)|) (b 2 ) = i?^ f (i? 2 ), where n is some non-zero power 



(or logarithmic for a purely Kolmogorov spectrum, applicable in the completely kinematic 
regime), applies only when the first two terms on the right of ([3]) approximately balance 
the next three, or equivalently, when the sum of the the first three terms on the right of 
(|i~0l) balance the first two terms on the right of (|TT|), since d t h = d t bW + o>*b«. But this 
state is very short lived because the small-scale field rapidly grows to near equipartition 



with the turbulent energy on a time scale ~ t c , and the first three terms on the right of (10) 
quickly dominate. These terms lead to an MHD turbulent spectrum with most of the energy 
contained on the outer scale of the turbulence. In principle, the mean field can still slowly 
grow, while the combination of field line stretching and nonlinear damping in a turbulent 
cascade maintain the small-scale field at a steady energy density. A key question ( [Field 1995]) 



is: what is the integrated kinetic vs. magnetic energy on the outer scale of the turbulence? 
We believe that this ratio is essential in assessing whether the dynamo works, and suggest 
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that a dynamo can work as long as there is a mismatch. 

Finally, note that the growth of the small-scale field necessarily means that even galaxies 
at redshifts > 2 would exhibit substantial fields when observed by Faraday rotation. This is 
because a typical turbulent energy-containing eddy in a galaxy is ~ lOOpc across, so that any 
line of sight through a disk will have significant RMS mean field component (Blackman 1998). 
A rough lower limit comes from assuming that the small-scale field builds up to equipartition 
with the turbulent motions (~ few x 10 _6 G in our Galaxy). Even when observed edge on, 
there are only ~ 100 such cells along the line of sight. The observed mean field would then 
be ~ 10 times smaller, or > 10~ 7 G. The mere existence of apparent mean fields is therefore 
not a strong test for dynamo theory even for large redshift Galaxies. A better observational 
test would be to measure the scale and pattern of field reversals at these redshifts. The 
existence of small-scale fields is guaranteed to accompany any large-scale field produced by 
dynamo action and does not preclude its operation. 
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